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Abstract 

This  paper  studies  a  general  class  of  two-dimensionaJ  systems  of  the  cubic  nonlineeir 
Schrodinger  type  (2DNLS),  defined  by  idtq  +  Oiq  =  pq  and  O2P  =  Oz{q*q),  where  each 
On  =  D^^^didj.,  n  =  1,2,3,  is  a  lineeir,  second-order,  operator  with  constant  coefficients. 
This  class  generahzes  the  Djordjevic-Redekopp  (DR)  system,  which  has  previously  been 
encountered  in  the  context  of  water  waves.  Integrability  is  characterized  siirply  in  terms 
of  covariant  conditions  on  the  0„.  We  obtain  all  int^grable  cases,  including  the  known 
cases  Davey-Stewartson  I  and  II  as  well  as  other  known  integrable  cases.  The  2DNLS  is 
modulationally  stable  if  >  0  Vfc,  where  =  kikjD^jK  All  other  regimes 

are  modulationally  \mstable  and  have  projections  satisfying  the  ordinary  (ID)  NLS  with 
soliton  solutions,  though  in  all  known  cases  these  ID  solitons  are  imstable  with  respect  to 
transverse  perturbations.  The  self-focusing  regime  is  characterized  by  the  eigenvalues  of 
the  d\J2:  Ox  and  O2  must  both  be  elliptic,  and  for  that  choice  of  variables  for  which 

and  both  have  positive  signature,  must  have  at  least  one  negative  eigenvalue.  The 
self-focusing  regime  is  distinct  from  the  modulationadly  stable  regime  and  also  from  the 
integrable  regime,  while  the  integrable  cases  may  be  modulationally  stable  or  unstable. 
There  are  no  soliton  solutions  known  in  those  integrable  cases  that  are  modulationally 
stable,  whereas  those  integrable  cases  in  which  2D  solitons  «u-e  known  correspond  to  the 
modulational  instability  regime. 
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1.  Motivation 


Reduction  of  nonlinear  physical  problems  to  systems  of  evolution  equations  with 
known  properties  has  proven  to  be  quite  fruitful  in  recent  years.  Though  a  comprehensive 
bibliography  is  beyond  the  scope  of  the  present  paper,  the  flavor  of  such  computations 
may  be  seen  in  our  o'wti  work  on  (l-t-l)-dimensional  “intense  propagation”  problems  as 
they  arise  in  plasma  physics  (see  Kates  ic  Kaup  1989a, b;  1991;  1992a;  1992b  and  refer¬ 
ences  cited  therein)  and  cosmology  (Kates  &  Kaup  1988).  Useful  examples  may  also  be 
found  in  a  series  of  papers  on  the  “reductive  perturbation  method”  (Ichikawa  et  al.,  1976; 
Taniuti,  1974).  Considering  the  wide  variety  of  applications,  even  modest  improvement 
in  our  imderstanding  of  the  resulting  nonlinear  evolution  equations  could  be  of  enormous 
value. 

If  nonlinear  efl'ects  are  not  too  strong,  and  if  s>Tnmetry  permits  the  assumption  of 
dependence  on  only  one  spatial  coordinate,  then  the  essential  features  of  a  nonlinear  pulse 
propagation  problem  may  often  be  deduced  from  the  well-known  properties  of  the  one¬ 
dimensional  nonlinear  Schrodinger  equation  (NLS) 


iqt  +  qxx±\q\^q~0  .  (1) 

In  the  reduction  of  a  complicated  ID  system  to  the  NLS,  the  most  important  computation 
is  the  dependence  of  the  sign  in  (1)  on  the  parameters  of  the  problem:  The  minus  sign 
indicates  modulational  stability,  while  the  plus  sign  indicates  modulational  instability  and 
soliton  solutions. 

In  this  paper,  we  investigate  the  qualitatively  much  richer  range  of  phenomena  asso¬ 
ciated  with  “2D”  instabilities  arising  in  nonlinear  pulse  propagation.  That  is,  we  consider 
the  evolution  of  nonlinear  instabilities  having  functional  dependence  on  (2-1-1)  dimensions 
(two  spatisd  dimensions  and  time). 

An  important  cirea  of  application  concerns  2D  modulations  of  electromagnetic  pulses. 
Since  Eq.  (1)  governs  ID  modulations  of  weakly  relativistic  EM  pulses,  one  might  have 
supposed  that  2D  modulations  of  EM  waves  could  be  described  simply  by  replEw:ing  the 
ID  dispersion  term  g**  in  the  NLS  (1)  with  a  2D  dispersion  term: 


iqt  -1-  [D^xdl  -f  2D^yd^dy  -f  Dyydl]q  ±  \qfq  =  0  ,  (2q) 


where 


D  — 


(26) 


A  system  of  the  form  (2)  is  indeed  obtained  for  a  circularly  polarized  beam  if  longitudi¬ 
nal  perturbations  are  negligible  (Spatschek,  1977).  This  case  corresponds  to  the  “thin- 
beam”  approximation.  However  in  many  applications  both  longitudinal  and  transverse 
perturbations  may  occur.  Careful  analysis  of  fluid  modes  associated  with  longitudineil  and 
transverse  ponderomotive  forces  (Kates  &  Kaup,  1993)  shows  that  a  proper  description  of 
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2D  modulations  of  EM  pulses  then  involves  two  independent,  coupled  potential  functions. 
This  point  was  recognized  by  Karpman  (1990). 

Some  hints  as  to  generic  2D  behavior  are  provided  by  water  waves  (Djordjevic  k. 
Redekopp,  1977  (DR  in  what  follows);  Ablowitz  &  Segur,  1979  (“AirS’’  in  what  follows)): 
CapiUary- gravity  waves  «ire  described  by  the  “DR”  system,  which  (for  historical  reasons) 
is  expressed  in  the  form 


+  ^r,r,  =  -PiA*A)^  .  (36) 

Both  (1)  and  (2)  are  special  cases  of  (3):  If  the  initial  data  depend  only  on  ^  or  r/  or 
some  lineeu'  combination  thereof,  the  DR  system  (3)  may  be  reduced  to  the  ID  nonlinear 
Schrodinger  equation  (1).  By  properly  choosing  the  coordinates,  one  may  write  (2)  in  the 
DR  form  (3),  in  which  case  =  0,  so  that  the  potentials  are  decoupled.  Eq.  (2)  or 
equivalently  (3a)  with  xi  =  0  will  be  refered  to  below  as  the  “decoupled  2DNLS.” 

As  written,  the  DR  system  (3)  contains  six  “arbitrary”  constants:  a,  0,  A  /x,  Xi 
Xi  •  (However,  by  suitable  coordinate  transformations,  one  can  reduce  the  dimensionality 
of  this  piirameter  space  to  two  parameters  and  two  signs,  as  is  occasionally  convenient.) 
Since  a  rather  large  class  of  2D  systems  with  cubic  nonlinearity  can  be  reduced  to  DR  form, 
it  represents  an  important  advance  in  generalizing  the  NLS  to  two  dimensions!  However, 
as  discussed  below,  the  DR  system  still  does  not  include  all  possible  2D  generalizations  of 
the  NLS.  Indeed,  it  does  not  contain  all  possible  2DNLS  integrable  cases.  Moreover,  in 
practice,  reduction  of  a  given  system  to  the  DR  form  (3)  may  require  a  tediou^  search  for 
appropriate  coordinate  transformations.  These  considerations  already  suffice  to  motivate 
our  introduction  of  a  more  general,  covariant  form  (Section  2),  which  seems  to  deserve  the 
designation  “c«inonical.”  This  canonical  form  includes  a  much  wider  class  of  2D  general¬ 
izations  of  the  NLS  and  also  facilitates  a  comprehensive  treatment  of  rJl  integrable  forms 
of  the  2DNLS.  The  importance  of  this  form  was  first  noted  by  Zakheirov  and  Schulman 
(1990). 

In  Section  3,  we  compute  the  regimes  of  modulational  instability  for  the  generic 
2DNLS.  More  precisely,  we  study  perturbations  about  spatially  constant,  bu|  nonlinear 
solutions.  Modulational  instability  of  the  2DNLS  is  always  associated  with  the  existence 
of  one-dimensional  NLS  soliton  solutions.  Whether  such  ID  solitons  truly  arise  in  a  system 
depending  on  two  spatial  coordinates  is  of  course  a  separate  question. 

The  canonical  form  turns  out  to  be  of  great  benefit  in  understanding  self-focusing. 
Self-focusing  has  been  studied  for  those  regions  of  the  DR  pareimeter  space  accessible  to 
water  wave  systems:  A&S  showed  that  a  singularity  in  moment  of  inertia  occurs  in  their 
“region  F,”  charcicterized  by  xi,  0-,  ct,  ^  positive,  x,  and  *'  =  X~Xii3/oi  negative  in 

Eqs.  (3).  They  also  studied  quasi-self-similar  solutions  in  this  region,  obtaining  growing 
solutions  consistent  with  a  self-focusing  interpretation.  The  question  of  self-focusing  in 
their  “Region  E”  (xi,  0,  Q,  /x,  A,  x,  and  v  all  positive)  was  left  open.  (Incidentally,  the 
A&S  paper  contains  a  few  minor  misprints:  The  quantity  A^  given  in  their  coordinate 
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transformations  should  read  |A1  throughout.  In  their  Eq.  (4. Id),  Xi  should  be  replaced  by 
Xi-) 

However,  several  intended  applications  to  experiments  in  plasma  physics  and  nonlinear 
optics  require  at  the  very  least  a  knowledge  of  self-focusing  conditions  for  the  entire  r«inge 
of  DR  parameters  (not  just  the  region  of  DR  parameter  space  accessible  to  water  waves). 
For  example,  the  interesting  question  of  self-focusing  for  positive  x  and  negative  v  did 
not  arise  in  the  A&S  paper,  because  for  water  waves  v  and  x  have  the  same  sign.  Other 
apphcations  may  require  self-focusing  results  for  systems  that  cannot  be  reduced  to  DR 
form  but  that  can  be  put  into  our  canonical  2DNLS  form:  for  example,  (Kates  &  Kaup, 
1992b)  the  Zakharov  equations  (Zakharov,  1972). 

In  Section  4,  we  derive  self-focusing  conditions  for  the  general  2DNLS  in  terms  of 
the  signatures  of  three  operators.  As  a  consequence  we  will  be  able  to  answer  the  open 
question  referred  to  by  AirS  concerning  the  possibility  of  self-focusing  for  their  “region  E” 
in  the  negative.  On  the  other  hand,  self-focusing  is  possible  for  positive  x  and  negative  u 
(eissuming  the  other  pcirameters  are  positive). 

2.  Canonical  Form  for  the  2DNLS  and  Integrability 


Of  pcirticiilar  interest  for  applications  such  cts  multi-dimensional  plasma  phenomena  is 
the  question  of  whether  par«imeters  can  be  chosen  such  that  the  system  is  integrable.  In  this 
case  the  system  may  have  nontrivial  soliton  solutions.  In  this  section,  we  discuss  a  general 
canonical  form  for  two-dimensional  generalizations  of  the  NLS  (“2DNLS”).  Conditions  for 
integrability  are  most  simply  and  conveniently  expressed  in  terms  of  this  canonical  form 
instead  of  the  DR  form. 

We  begin  our  discussion  by  observing  that  if  both  eigenvalues  of  Dij  have  the  same 
sign  as  the  nonlinear  term,  then  we  may  rescale  the  decoupled  NLS  (2)  to 

idtq  -I-  -i-  \q\^q  =  0  .  (4) 

In  two  dimensions,  Eq.  (4)  is  of  considerable  interest  in  its  own  right,  because  it  exhibits 
self-focusing  behavior  (Synakh,  1975;  Spatschek,  1977).  Eq.  (4)  would  appear  to  be  the 
“simplest”  possible  multi-dimensional  generalization  of  (1)  in  any  number  of  dimensions; 
if  the  ID  version  of  a  physical  system  obeys  the  NLS  (1),  then  it  sometimes  happens  that 
the  multi-dimensional  version  can  be  reduced  to  (4).  However,  its  properties  are  quite 
different:  In  more  than  one  dimension,  Eq.  (4)  is  not  integrable. 

All  known  integrable  forms  of  the  multidimensional  NLS  involve  an  auxiliary  function. 
Such  a  function  arises  naturally  in  the  physics  of  water  waves  (DR),  in  the  Zakharov  system 
(see  below)  and  in  electromagnetic  pulses  propagation  (Kates  and  Kaup,  1992b). 

Recently,  Boiti,  Pimpinelh  and  Sabatier  (1992)  demonstrated  that  there  are  only  three 
distinct  integrable  forms  of  a  2DNLS.  Below,  we  will  show  that  all  three  may  be  expressed 
compcictly  in  terms  of  the  “canonical”  form  (Zakharov  and  Schulman,  1990) 
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idtq  +  Oiq  =  pq  (5a) 

O2P  =  03(9*9)  ,  (55) 

where  each  of  the  operators  0„  in  (5)  is  a  second-order  linear,  symmetric  (dispersion) 
operator  of  the  form 


On  =  Di^Jdl  +  2Dil^d,dy-\-Dl;Jdl  .  (5c) 

The  criteria  for  integrability  of  (5)  now  become  simply  that  at  least  one  of  O2  and  O3  not 

be  elliptic, 

£>;■'  =  KD\f  (6.1) 

be  satisfied  for  some  constcint  k, 

=  0  ,  (66) 


and 

det(D(”>)  #  0  (6c) 

for  each  of  n  =  1,2,3.  Thus,  expressing  the  system  in  the  canonical  form  (5)  offers  the 
advantage  of  allowing  a  straightforward,  coordinate-independent  test  of  integrability. 

Note  that  whenever  det(j9^”))  =  0,  then  is  simply  a  projection.  Thus  condition 
(6c)  ensures  that  each  0„  remains  two-dimensional.  Condition  (6b)  means  that  in  some 
sense,  the  operator  O2  must  be  orthogonal  to  O3  (and  naturally  also  to  Oj). 

We  now  relate  conditions  (6)  to  the  integrability  conditions  of  Boiti,  Pimpinelli  and 
Sabatier  (1992);  Using  powerful  an«ilytical  techniques,  they  showed  that,  in  multiple  di¬ 
mensions,  the  only  possible  integrable  equations  having  at  most  two  derivatives  are  ob¬ 
tainable  from  a  Lax  peiir  of  the  following  form: 


d2Vi  =  qv2  (7a) 

diV2=rvi  (76) 

idtVi  =  aidlvi  -H ibidivi  +  +  eiV2  (8a) 

idtV2  =  U2dlv2  -F  162^^2  -b  C2diVi  +  d2Vi  +  e2V2  ,  (86) 


where  di  =  d/dwi;d2  =  d/dw2‘,  u^i  and  W2  are  independent  spatial  variables;  Uj  and  V2 
are  eigenfunctions;  ai  and  <T2  *ire  constants;  and  6i,62,ci,C2,di,d2.ei  and  62  are  functions 
to  be  determined.  Eqs.  (7)  constitute  a  scattering  problem  which  was  first  studied  by 
Ablowitz  and  Habermcm  (1975).  For  given  q  and  r,  particular  scattering  data  are  deter¬ 
mined.  Conversely,  given  the  scattering  data,  one  can  reconstruct  the  potentials,  9  and 
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The  Lax  pair  (7)  «ind  (8)  can  only  have  a  common  solution  when  certain  integrability 
conditions  are  satisfied.  These  integrability  conditions  yield  precisely  the  integrable  forms 
of  the  2DNLS  detailed  by  Boiti,  Pimpinelli  and  Sabatier  (1992).  One  obtains  them  by 
cross-differentiating  (7)  and  (8),  using  (7)  and  (8)  to  simplify  the  resulting  conditions. 
One  finds 


Cl  =  029  (9a) 

C2  ■=  a iT  (9b) 

^6i  =  0  (9c) 

dibs  =  0  (9d) 

Cl  =  iqi>2  -  <^2^29  (9e) 

ds  =  trbi  —  Oidir  (9/) 

and  the  evolution  equations 

i{dt  -  bidi  -  b2d2)9  =  (<^idi  -  a^dDq  -f  (di  -  cs  -I-  td2b2)9  (10a) 


-i(9t  -  bidi  -  b2d2)r  =  (oidj  -  osdDr -I- (di  -  es  -  *dibi)r  ,  (10b) 

where  the  potentials  di  and  cs  are  detennined  by 

dsdi  =  — 2oidi(r9)  (Ha) 

dies  =  -2o2d2(r9)  .  (11b) 

The  quantities  bi  and  bs  are  trivial,  being  simply  a  group  velocity  for  the  q  and  r  fields. 
Thus  we  take  them  to  be  zero.  Now,  if  we  identify 

r  =  <To9*  ,  (12) 

where  oq  ~  ±1,  then  the  two  equations  (10)  are  obviously  equivalent.  Defining 

p  =  di  -  62  ,  (13) 

we  then  have  from  (11) 

did^p  =  -2ao(oidf  -  02dl){q*q)  (14a) 

and  (10a)  becomes 

idtq  =  (<71^1  -  <^2^2)9  +  P9  •  (14b) 
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Lastly,  if  we  define 


02  =8182 

(15a) 

O3  =  -2aQ{(Ti8i  -  0282) 

(156) 

Oi  =  (cT^af  -<72d|)  , 

(15c) 

the  system  takes  the  canonical  form  (5)  and  satisfies  the  conditions  given  in  (6).  Conversely, 
conditions  (6)  are  invciriaint  under  general  lineair  coordinate  transformations  as  well  as  imder 
constant  rescalings  of  the  dependent  variables.  This  invariance  simplifies  considerably  the 
problem  of  determining  when  an  equation  of  the  form  (5)  is  integrable.  One  does  not  have 
to  secirch  for  special  coordinate  scalings.  One  simply  checks  whether  or  not  conditions 
(6)  hold  in  some  coordinate  system.  Wlien  (6)  are  satisfied,  then  one  can  always  find 
coordinates  such  that  the  0„  are  given  by  (15). 

Let  us  now  make  some  remarks  on  this  canonical  form.  First,  there  is  an  arbitrariness 
in  the  choice  of  signs  of  wi  and/or  W2-  By  using  this,  we  could  absorb  the  sign  oq  into 
either  toi  or  ui2-  Note  that  (7)  would  be  transformed  if  the  sign  of  oq  is  so  absorbed. 
However,  it  is  convenient  to  carry  oq  in  the  computation  and  leave  the  coordinates  fixed. 
Furthermore,  the  dispersion  operator  O2  is  not  necessarily  restricted  to  be  hyperbolic,  as 
it  may  appear  at  first  glance.  When  the  dispersion  operator  O2  is  hyperbolic,  then  it  is 
only  necessary  to  choose  the  coordinates  wi  and  W2  to  be  real.  Then  in  order  to  obtain 
real  operators  O3  and  Oi ,  one  would  choose  ai  and  02  real.  In  this  case,  O3  and  0%  may 
be  either  elliptic  or  hyperbolic,  depending  on  the  signs  chosen  for  aj  and  02. 

However,  suppose  the  actuzil  problem  demands  O2  to  be  elliptic.  One  may  also  treat 
this  case  by  allowing  u;i  and  W2  to  be  complex.  (In  the  above  derivation,  di  and  82 
were  only  required  to  be  linear  operators.  They  quite  easily  could  have  been  a  linear 
combination  of  derivatives  of  read  coordinates  with  complex  coefficients  as  below.)  The 
simplest  example  is  to  take 


di=dx+idy  (16a) 

d2=dx-idy  ,  (166) 

where  x  aind  y  axe  reail  coordinates.  Now  the  dispersion  operators  all  have  to  be  real  for 
physical  reasons,  so  care  has  to  be  taiken  to  ensure  that  this  remains  so  when  and  W2 
are  complex.  The  condition  that  O3  (and  Oi)  be  real  requires  02  =  —<^1,  so  that  only  one 
of  <Ti  and  <72  may  be  chosen  arbitrarily.  One  may  show  that  for  O2  elliptic,  O3  (and  Oi) 
must  be  hyperbolic.  This  is  in  contrast  to  the  case  O2  hyperbobc,  in  which  O3  (and  Oi ) 
could  be  either  elliptic  or  hyperbolic.  Also,  one  may  verify  that  this  does  not  affect  the 
choice  (12),  which  remains  a  valid  identification. 

This  case  of  O2  elliptic  contains  the  case  of  the  classical  Davey-Stewartson  equation 
(which  has  been  called  the  “Davey-Stewartson  11” )  in  hydrodynamics.  We  refer  the  reader 
to  the  literature  [Kaup  (1980),  Fokas  and  Ablowitz  (1983),  Fokas  and  Ablowitz  (1984), 
Boiti,  Leon,  Martina  and  Pempinelli  (1988),  Boiti,  Leon  and  Pempinelli  (1989),  Boiti, 
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Pempinelli  and  Sabatier  (1992),  Arl-idiev,  Pogrebkov  and  Poliv-anov  (1989)]  for  the  solution 
of  the  initial  value  problems  and  various  specific  solutions. 

Now,  an  important  point  in  regard  to  this  canonical  system  is  that  any  system  which 
can  be  expressed  in  the  DR  form  (3)  can  also  be  expressed  in  the  canonical  form  (5).  (The 
converse  is  not  always  true.)  The  transformation  can  be  easily  obtained  by  differentiating 
(3b)  with  respect  to  (  and  making  the  identifications 


O2  =  — -\ - 

(17a) 

Xi  Xi  ^ 

(176) 

Oi  =  -1-  A9| 

(17c) 

q  =  A 

(17d) 

p  =  x{A' A)  +  X\h^  ■ 

(17e) 

(Note  that  x  corresponds  to  tj  and  y  to  The  “decoupled”  case  of  vanishing  Xi  is  excluded.) 

The  so-called  Davey-Stewartson  I  equations  (DS-I) 

-iQr=\{Q..+Qyy)-^Q{4',±Q^Q)  (18a) 

{<l>r.-<t’yy)±nQQ‘)y=0  (186) 

and  DS-II  equations  (for  the  upper  sign) 

-iQr  =  \iQ^.  -  Qyy)  '  Q(0y  ±  Q* Q)  (19a) 

{(f'xx  +  <t>yy)  i  AQQ*)y  =  0  (196) 

are  the  integrable  cases  of  DR. 

Conditions  for  integrability  of  a  system  expressed  in  the  DR  form  (3)  may  be  obtained 
by  the  foUowing  tedious  procedure:  one  introduces  constant  transformations  to  normal¬ 
ize  all  but  two  of  the  coefficients  in  (3).  Inspection  of  the  resulting  equations  reveals 
that  reduction  to  Eqs.  (18)  (DS-I)  requires  positive  A  and  satisfaction  of  two  additional 
conditions: 

-  -2xA  (20a) 

and 

Qp  =  -A  .  (206) 

If  A  is  negative,  the  DR  system  (3)  reduces  to  Eqs.  (19)  (DS-II)  if  x  is  positive  and  (20) 
hold.  Applying  (6)  to  the  operators  On  defined  in  Eqs.  (17a-c),  one  sees  that  the  conditions 
for  integrability  are  consistent  with  those  given  in  (20)  for  the  DR  form. 
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The  “standard”  forms  (18)  and  (19)  for  DS-I  and  DS-II  appear  to  distinguish  a  direc¬ 
tion  of  propagation  (longitudinal  coordinate)  from  a  transverse  coordinate,  whereas  there 
is  no  such  distinction  in  (5).  Now,  of  course  such  a  distinction  may  be  defined  in  the 
underlying  problem,  but  it  is  not  inherent  in  the  DS  equations  themselves,  since  one  can 
easily  transform  either  DS  system  to  (14). 

In  some  cases,  we  can  reverse  the  procedure  and  convert  (14)  to  DR  form.  However, 
we  emphasize  that  not  all  integrable  cases  can  be  expressed  in  DR  form:  Suppose  that  a 
transformation  exists  such  that  p  can  be  decomposed  into 

P  =  (”1^1 +«252)4>  ,  (21a) 

where  x  is  a  coupling  constant  and  {nidi  +  n2d2)  is  the  gradient  in  some  “transverse” 
direction.  From  (14a),  one  then  must  obtain  the  equation  for  $  in  the  form 

8182^  =  2'r{nidi +n2d2){q*q)  ,  (216) 

where  7  is  simply  a  sign.  One  finds  that  this  can  be  done  whenever  the  following  conditions 
can  be  satisfied: 


Cl  =  -'rni^  (22q) 

,  (226) 

where  x  will  then  be  given  by 

X  =  -47nin2  .  (23) 

When  O2  is  elliptic,  then  since  02  =  — tri*,  a  solution  will  always  exist  for  the  Ui's 
where  n2  =  nj*,  leading  to  a  real  transverse  gradient  operator,  {nidi  +  n282)  (recall  here 
that  the  p«irtials  are  complex  operators).  However,  in  the  case  where  O2  is  hyperbolic  and 
the  coordinates  wi  and  W2  are  re«il,  from  (22)  it  is  clear  that  a  solution  exists  only  if  O3  is 
elliptic  (DS-I).  When  in  addition  O3  is  hyperbolic,  then  no  real  solution  can  exist  for 
and  n2,  i.e.,  at  least  one  of  the  direction  cosines  would  have  to  be  complex.  Consequently 
the  DR  form  is  not  a  canonical  form  covering  all  possible  integrable  cases.  In  particular, 
it  cannot  cover  the  case  where  both  O2  and  O3  are  hyperbolic.  For  this  reason,  the  form 
(5)  is  preferable,  because  it  does  cover  all  integrable  cases. 

3.  Modulational  Instability  of  the  2DNLS 

As  in  the  one-dimensional  NLS  (1),  important  hints  to  the  asymptotic  behavior  of 
solutions  of  the  generic  2DNLS  (5)  can  be  obtained  by  studying  modulational  stability 
foUowing  the  method  of  Paper  III.  We  first  observe  that  (5)  admits  spatially  constemt 
solutions  of  the  form  p  =  po  =  const,  q  =  qo{t),  where  po  and  go  are  any  solutions  of 


idtqo  =  Poqo 


(24) 
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for  example, 


go  =  exp{ipoO 


(25) 


We  next  write 

P  =  Po(l+Pi)  (26a) 

9  =  9o(l  +  9i)  (266) 

(where  qi  and  pi  depend  on  all  variables)  and  linearize  Eqs.  (5)  in  pi  and  gi ,  obtaining 


idtqi  +  Oiqi  -poPi  =  0 

-idtq*  +  Oiqi  -poPi  =  0 

O2P]  =  +  9o) 

Po 


(27a) 

(276) 

(27c) 


Expressing  the  dependent  variables  qi,  q*,  and  pi  as  Fourier  integrals  and  being  careful 
to  include  the  possibility  that  the  time  dependence  of  the  phase  may  be  complex,  we  then 
obtain  the  dispersion  relation 


(28) 


where  is  defined  by 


exp(t^  •  x)  =  —On  txp{ik  ■  x)  ,  (29a ) 

that  is, 

.  (296) 

Let  us  assume  that  is  nonvanishing.  Since  it  can  be  shown  that  higher  wavenumbers 
are  stabilizing,  we  obtaun  the  condition 

p(i)p(2)p(3)  >  0  Vit  (30) 

for  modulational  stability  of  the  general  canonical  system  (5). 

One  important  application  of  (30)  concerns  one-dimensional  subspaces  (see  also  Dis¬ 
cussion):  For  any  fixed  k,  we  can  define  y  =  k-x  and  substitute  p  =  p(t,  y)  and  q  =  q{t,  y) 
into  (5).  Assuming  nonvanishing  (5b)  can  easily  be  integrated  twice,  and  p  can  be 
eliminated  form  (5a),  resulting  in  a  different  one-dimensional  NLS  (1)  for  each  k.  It  is  then 
easy  to  verify  that  (30)  yields  the  correct  NLS  modulational  stability  condition.  Additional 
consequences  will  be  discussed  in  the  final  section. 

It  is  especially  interesting  to  combine  (30)  with  the  integrability  conditions  (6):  these 
imply 

,  (31) 
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and  thus 


kP'2*  >  0  (32) 

is  the  condition  for  modulationail  stability  of  cin  integrabk  2DNLS  system.  Consequences 
will  be  discussed  in  the  final  section. 

Although  the  above  derivation  already  covers  the  DR  form  (3)  in  principle  (via  Eqs. 
(17)),  for  future  reference  we  repeat  the  argument  explicitly  for  the  DR  form  (3):  The 
spatially  constant  solution  for  A  is 

Ao(t)  =  aoexp(-txTao)  (33) 


with  $0  =  0.  We  take 


a  =  Ao(l +ai(T,C,Tj))  , 


and  hnearize  (3)  in  ui  and  The  dispersion  relation  is  then 


(34) 


u 


k^{Xs'm‘^{6)  +  /iCos^(^))[(Asin^(0)  +  ficos^{9))k^  -I-  2aoX 


2al3xi  sin^(g) 
cos'^{6)  +  osin^(0) 


(35) 

where  =  fccos(0),  k^  =  fcsin(0),  i.e.,  +  k^.  Since  evidently  higher  wavenumbers 

k  are  stabilizing,  the  maximal  instability  region  may  be  determined  by  considering  the 
limit  A:  0  of  7^/fc^  .  The  resulting  inequality  reduces  properly  in  the  limiting  cases  6  =  0 
and  6  =  7r/2  to  the  known  conditions  for  modulational  instability  of  the  “longitudinal’' 
and  “transverse”  NLS.  It  is  also  is  of  great  utility  in  understanding  the  changes  in  the 
qualitative  behavior  of  solutions  of  the  DR  system  for  limiting  values  of  the  coefficients. 


Finally,  we  note  that  the  zeroth-order  solution  (33),  which  was  assumed  spatially 
constant,  may  be  generalized  to  an  arbitrary  plane  wave  by  Galilean  invarieince  without 
changing  the  regimes  of  modulational  stability. 


4.  Self-focusing 


Probably  the  most  dramatic  difference  between  ID  and  2D  modulations  is  the  possi¬ 
bility  of  a  violent  “self-focusing”  singularity  (Zakharov  k  Synakh,  1975;  Spatschek,  1977). 
Self-focusing  is  also  a  well-known  phenomenon  in  nonlinear  optics  (Kelley,  1965).  For  the 
general  2DNLS  (5),  which  contains  two  independent  potentials,  the  situation  is  a  quite  a 
bit  more  complicated  than  for  the  decoupled  2DNLS  (4),  which  involves  only  one  potential. 
Below,  we  derive  restrictions  on  the  operators  On  in  (5)  and  give  conditions  on  the  initial 
conditions  leading  to  self-focusing. 

We  are  interested  in  the  time  evolution  of  the  generalized  “moment  of  inertia”  integral 


j=  f  [ 

J  J  .. 


XiXjq*q 


(35) 
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We  assume  that  Oi  and  Oi  are  elliptic  and  that  all  fields  vanish  at  infinity  sufficiently 
rapidly  in  what  follows  (see  end  of  this  section  for  discussion).  Without  further  loss  of 
generality,  we  take  the  matrix  to  have  two  positive  eigenvalues.  Note  that  in  the 
absence  of  singularities  J  is  a  positive  definite  quantity.  (This  would  not  be  true  for 
hyperbolic  Oi-)  For  our  purposes,  the  question  of  self-focusing  amounts  to  asking  under 
what  circumstances  J  develops  a  singularity  in  finite  time. 

It  may  be  verified  (see  Appendix  A)  under  the  above  zissumptions  that 


J 


=  8/i  , 


where 


J  J  d^df)  Hx 


(37a) 


Hi  =  +  \pq-q 


Consider  a  ( “generalized  DR” )  system  of  the  form 

iAr  -I-  OiA  =  xA*A^  -h 


(376) 


If  Ii  exists,  is  conserved,  cind  is  initi«illy  negative,  then  integration  of  (I  6)  implies 
that  J  will  develop  a  singularity  in  finite  time.  Now,  for  arbitrary  0„,  h  is  not  generally 
conserved,  even  for  diagonal  Oi  and  O2  (i.e.,  the  DR  system).  Nevertheless,  &s  we  shall 
now  show,  Ii  will  be  constant  in  time  imder  the  assumptions  stated  above. 


(38a) 


02^  =  -I3{A*A)^  , 


(386) 


where  Oi  and  O2  are  as  in  (5).  The  DR  form  (3)  corresponds  to  the  special  case  of  (38)  in 
which  Oi  and  O2  are  both  diagonal.  Our  procedure  will  now  be  to  write  the  2DNLS  (5) 
in  the  form  (38)  and  obtain  conservation  laws  for  /j  in  this  form. 


Let  us  write 


q  =  A 


Eq.  (5b)  becomes 


p  =  xA*A  -h  P 


O2P  =  (O3  -  x02)(AM)  . 


Since  O2  is  elliptic,  it  can  be  shown  (see  following  paragraph)  that  there  exists  a  linear 
transformation  to  coordinates  77,  ^  such  that 


O3  -  XO2  =  -bdl 
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for  some  constant  b,  where 


(43) 


3 

and  where  without  loss  of  generality  n  is  a  unit  vector.  Identifying 


(44) 


we  can  then  integrate  (41). 

Now  let  us  prove  the  above  statement  leading  to  (42):  Since  by  assumption  both 
eigenvalues  of  are  positive,  we  first  transform  it  into  the  identity  matrix  by  means 
of  a  consteint  linear  coordinate  transformation.  In  the  coordinate  system  (xi,X2),  we  thus 
seek  a  solution  of  the  matrix  equation 

■Olf  -  X^ij  +  =  0  .  (45) 

Eq.  (45)  is  equivalent  to  three  scalar  equations.  One  of  these,  the  trace  of  (45),  yields  an 
expression  for  b: 

i.=  -(DW  +  Dg*)  +  2x  .  (46) 

Let  rij  =  cos  6  and  n2  =  sin  6  for  some  6.  The  remaining  two  scalar  equations  then  yield 

X  =  5  (bI!*  +  d'S  ±  [(fll?’  -  Dg')’  +  4(r>'|')=]'/=‘)  (47) 


and 


cos 


20  = 


-  X 


D[\'  +  d'^)  -  2x 


<  1 


(48) 


Therefore,  imder  the  stated  assumptions  we  may  always  write  the  2DNLS  (5)  in  the  form 
(38). 

In  order  to  prove  conservation  of  I\,  it  is  convenient  to  define  a  Lagrangi2m  density 
for  (38): 


where 


L=-{ArA^-A*,A)-H{A,A\ii>) 


(49) 


"  =  E  +  §  E  +  XiAA'ii  .  (50) 

ij  ^  ij 

Equations  (38)  are  then  equivalent  to 


d^drj  L 


=  0 


(51) 


where  independent  Veu-iations  cire  to  be  performed  with  respect  to  A*,  A,  and  In 
particular,  (51)  implies 


d  /q 
dr 


=  0  , 


(52a) 
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where 


(526) 


where 


h  =  J  J  did,,  Ha  , 

(546) 

Ha=Yi  +  \(AA‘f  +  . 

(54c) 

On  the  other  hand,  substituting  (40)  and  (44)  into  (37)  and  taking  the  special  case  Xi  =  li 
we  find  that  Hi  =  H3.  Therefore,  Ii  is  conserved. 

Let  us  now  study  conditions  for  negative  Ii:  Since  the  first  term  in  Hi  is  positive 


definite,  self-focusing  requires  ; 

/4  <0  , 

(55a) 

where 

74  =  y  y  <fx  pq*q. 

(556) 

A  sufficient  condition  for  (55a)  is  that  O2  and  O3  have  opposite  signatures:  We  first  express 
p  and  q*q  as  2D  Fourier  integrcils 

P=  j  jd^k 

(56a) 

9*9  =  y  y  d'^k  pj^u^  , 

(566) 

where 

Uj- =  —  exp(ik  -  x)  . 

(57) 

We  then  have 

A  =  y  /  i’lPk  • 
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Observing  that  (5b)  is  a  linear  equation  in  p  and  q*q,  we  can  now  solve  for  V’jf  in  terms 
of  pjf  using 


•  (59) 

ij  O 

The  integrand  in  (58)  is  therefore  negative  definite  if  the  signatures  of  O2  and  O3  are 
opposite,  and  therefore  (55a)  is  satisfied  as  claimed. 

For  self-focusing,  we  are  interested  in  initial  data  such  that  not  just  1^,  but  also  Ii  is 
negative.  It  is  intuitively  clear  that  this  can  always  be  accomplished  if  q  is  chosen  to  have 
sufficiently  large  power  at  small  k. 

We  observe  that  if  the  signatures  of  O2  and  O3  are  the  same,  then  I4  is  positive 
definite.  In  this  case,  the  system  is  clearly  outside  of  the  self-focusing  regime  as  defined 
by  the  prediction  of  a  singularity  in  J.  This  is  the  case  of  “Region  E”  defined  by  A&S  for 
the  DR  system  (3). 

In  the  case  of  I?-  having  one  negative  and  one  positive  eigenvalue  (hyperbolic  O3), 
Eq.  (59)  no  longer  guarantees  opposite  signs  of  V’jE  and  pj  for  arbitrary  k.  An  example 
is  the  case  of  positive  x  negative  v  in  (3)  (other  parameters  positive)  mentioned  in 
the  introduction.  Nevertheless,  self-focusing  can  still  occur  if  the  initial  conditions  have 
enough  power  at  low  wave  number  and  also  most  of  the  power  in  q*q  is  concentrated  in 
modes  with  wave  numbers  k  intersecting  a  hyperbola  D\fkikj  =  -1  (i.e.,  in  the  notation 
of  (29a),  modes  satisfying  <  0). 

At  the  beginning  of  this  section,  we  assumed  that  Oi  and  O2  are  elliptic.  Now,  the 
condition  Oi  elliptic  obviously  is  needed  to  argue  for  a  singuleuity  in  J  (Eq.  (35)).  The 
condition  that  O2  must  be  elliptic  was  derived  for  the  DR  system  (3)  by  A&S,  and  their 
arguments  are  simply  repeated  here  for  the  2DNLS  (5);  First,  suppose  one  has  initial  data 
for  q  which  decays  sufficiently  rapidly  at  infinity  or  has  compact  support.  Then  partial 
integration  of  global  integrals  involving  only  q  and  its  derivatives  leave  no  boundary  terms. 
If  O2  is  elliptic,  then  appropriate  boimdary  conditions  for  (5b)  are  that  p  vanish  at  infinity, 
so  integration  by  parts  of  global  integrals  involving  p  also  leaves  no  boundary  terms.  This 
argument  breaks  down  if  O2  is  hyperbolic,  since  boundary  terms  involving  p  do  not  in 
general  vanish  at  infinity. 

Summzirizing,  for  the  self-focusing  regime,  the  operators  O2  and  Oi  in  (6)  must  be 
elliptic,  while  must  have  at  least  one  negative  eigenvalue  (taking  and  positive 
definite). 


5.  Discussion 

The  canonical  form  (5)  of  the  2DNLS  has  proven  to  be  quite  useful  in  understanding 
integrability,  modulational  instability  and  self-focusing.  An  overview  is  given  in  Table  1. 
[For  example,  the  entry  “DSI”  under  “Integrability”  in  the  sixth  row  of  the  table  means 
that  for  (Oi,  O2,  O3)  =  (elliptic,  hyperbolic,  elliptic),  conditions  (6)  are  equivalent  to  DSI.] 


15 


We  have  obtained  all  of  the  integrable  cases  (see  conditions  (6)  which  also  contain  some 
degenerate  cases  [Zcikharov,  1980]).  Eq.  (32)  is  the  condition  for  modulational  stability 
of  an  integrable  system.  Observe  that  the  DS-I  system,  which  is  known  to  admit  soUtons, 
is  modulationally  unstable,  as  one  would  have  expected,  while  the  DS-II  system,  which  is 
not  known  to  admit  solitons,  is  in  fact  modulationally  stable  (for  k  >  0).  It  thus  seems 
plausible  that  all  integrable  cases  that  do  not  give  rise  to  solitons  are  in  the  modulationally 
stable  regime. 

Reccill  that,  for  the  self- focusing  regime  as  we  understand  it  here,  the  operators  O2  and 
Oi  in  (6)  must  be  elliptic,  while  must  have  at  least  one  negative  eigenvalue  (taking 
jD(2)  £)(,!)  positive  signature  without  loss  of  generahty).  Hence,  the  regimes 

of  self-focusing  and  modulational  stability  (30)  are  mutually  exclusive.  For  applications  to 
practical  experiments,  it  is  especially  interesting  to  compare  (28),  the  dispersion  relation 
for  modulational  instability  (arbitrary  k),  with  the  conditions  for  self-focusing  (negativity 
of  /i;  see  (37))  in  light  of  Eq.  (59):  In  (28),  higher  wavenumbers  are  stabilizing,  and  in  /j, 
too  much  power  at  high  wavenumbers  prevents  self-focusing,  even  if  I 4,  (see  Eqs.  (55)  and 
(58))  is  negative. 

In  psu-ticulaj,  for  hyperbolic  O3,  self-focusing  requires  not  only  sufficient  power  at  low 
wavenumbers,  but  also  that  most  of  the  power  in  g*q  be  concentrated  in  modes  k  such 
that  <  0  (in  the  notation  of  (29a)).  These  are  just  the  unstable  modes  according  to 
(28). 

Next,  we  observe  that  the  conditions  for  integrability  and  those  for  self-focusing  are 
mutually  exclusive.  Suppose  the  contrary:  As  stated  above  Eq.  (6a),  for  integrability  at 
least  one  of  O2  and  O3  must  not  be  elliptic.  For  self-focusing,  O2  must  be  elliptic,  so  O3 
would  have  to  be  hyperbolic.  But  then  by  (6a)  Oi  would  also  be  hyperbolic,  which  violates 
our  conditions  for  self- focusing. 

In  Section  3,  we  saw  that  for  any  ID  subspace  defined  by  y  =  it  •  x  (where  k  is  any 
direction  such  that  does  not  vanish),  the  2DNLS  (5)  may  be  reduced  to  Eq.  (1), 
the  (ID)  NLS.  Modulational  instabilities  of  the  2DNLS  predicted  by  (30)  are  therefore 
associated  with  ordinary  NLS  solitons.  However,  from  the  work  of  A&S  we  know  at  least 
for  the  DR  equations  that  all  ID  solitons  are  unstable  with  respect  to  long-wavelength 
perturbations  in  the  remaining  spatial  coordinate  if  these  perturbations  are  compatible 
with  the  boundciry  conditions.  An  intriguing  question  for  future  work  is  thus  whether 
all  ID  .olitons  are  unstable  for  the  full  2DNLS  (5)  and,  if  so,  what  they  evolve  into. 
Of  special  interest  is  the  regime  where  all  three  operators  are  elliptic  and  has  two 
negative  eigenvalues  (taking  and  positive  definite).  This  regime  corresponds  to 
“Region  F”  of  A&S  in  the  DR  case  and  is  the  most  favorable  for  self-focusing.  Begging 
the  reader’s  tolerance  for  pure  speculation,  we  suggest  that  if  “nezuly  ID”  solitons  do  form 
and  persist  for  some  time  in  this  regime,  they  may  decay  due  to  growth  of  long- wavelength 
“transverse”  perturbations  and  eventucdly  develop  self-focusing  singularities. 

It  is  also  evident  from  (30)  that  a  system  which  reduces  to  a  modulationally  stable 
NLS  in  a  ID  subspace  defined  hy  y  =  k  •  x  may  still  be  imstable  in  2D.  In  particular,  it 
may  self-focus  or,  if  it  happens  to  b.-:  integrable,  form  2D  solitons.  This  simple  observation. 
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together  with  the  instability  of  ID  solitons  discussed  in  the  preceeding  paragraph,  poses  a 
severe  limitation  on  the  validity  of  conclusions  drawn  from  a  one-dimensional  analysis  of 
any  intrinsically  multi-dimensional,  NLS-type  system. 


Table  1:  Self-Focussing,  Modulational  Stability, 
and  Integrability  of  2DNLS  (Eqs.  (5)) 


Ox 

O2  O3 

Self¬ 

focussing 

Modulational 

stability 

Integra¬ 

bility 

H 

H 

H 

7 

U 

YES/NO 

H 

H 

E 

7 

U/S 

NO 

H 

E 

H 

7 

U/S 

DSIl/DSir/NO 

H 

E 

E 

7 

u 

NO 

E 

H 

H 

7 

U/S 

NO 

E 

H 

E 

7 

u 

DSI/NO 

E 

E 

H 

YES 

u 

NO 

E 

E 

E 

YES/NO 

U/S 

NO 

H  =  hyperbolic,  E=elliptic,  U  =  imstable,  S  =  stable 

Where  there  is  a  choice,  which  case  will  occur  wiU  depend  on  the  exact  values  of  the 
coefficients  (see  text). 

DSI,  DSII,  etc.:  (see  text) 
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Appendix  A:  Derivation  of  £q.  (42) 
Consider  the  moment  of  inertia  integral 


J  = 


(>11) 


We  assume  that  0\  and  O2  are  elliptic,  all  fields  vanish  at  infinity  sufficiently  rapidly,  and 
that  q*q  and  p  can  be  expressed  as  2D  Fourier  integrals. 

Eq.  (5a)  may  be  expressed  in  the  form 


idtq  = 

{A2a)  , 

where 

"H  =  +  p 

{A2b) 

is  Hermiticin.  Using  (A2),  one  obtains 

dt  qxVq]  . 

(A3) 

Differentiating,  one  obtains 

Jd^x(-q-[n  ,f-V],  +  c.c.)  , 

(A4) 

where  [  ,  ]  indicates  the  commutator  of  two  operators.  We  note  the  commutation  relations 

[-O1  ,x.V]  =  -20i 

(A5) 

and 

[p  ,x  •  V]  =  —X  •  V 

(A6) 

Peutially  integrating,  one  obteiins 

^=8//d»x^DljV-4/,  . 

{Ala) 

where 

h  =  J  J  (Pxq*q{x  •  Vp)  . 

{Alb) 

Let  us  write  p,  q*q,  and  x  •  Vp  as  2D  Fourier  integrals 

P  =  J 

(A8) 
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where 


9*9  =  y  j  PkH  - 


«jt  =  ^  exp(tfc  •  x) 


/s  then  reduces  to 


h  =  JJ  dH  . 


(AlO) 

(All) 

(A12) 


For  simplicity,  we  drop  the  subscripts  on  u,  rp,  and  p  in  what  follows.  Fi-om  (A9) 
-  (AlO)  one  obtains 


(A13) 


Now,  since  O2  is  elliptic,  Eq.  (5b)  implies 


Ip  =  Rp  , 

(A14a) 

where 

p(3) 

^  -  x>(2) 

(A146) 

and 

(A14c) 

Using  (A13)  and  (A14)  and  vigorously  integrating  by  parts  in  fc-space. 

one  obtains 

h  =  -J  jd^k  tptpj^  . 

(A15) 

Inverting  the  Fourier  transforms,  one  then  has 

h  =  -  j  J  dFx  pq*q  . 

(A16) 

Combining  (A16)  and  (A7),  one  then  obtains  Eq.  (42). 

Note  that  ellipticity  of  O2  was  essential  to  the  argument.  Otherwise,  the  quantity  R 
defined  in  (A14b)  wotild  have  been  singular. 

Acknowledgements 

This  research  was  supported  in  part  by  the  Office  of  Naval  Research  and  the  Air  Force 
Office  of  Scientific  Research. 


REFERENCES 


Ablowitz,  M.J.  &  Haberman,  R.  1975  Phys.  Rev.  Letters  Zh,  1185. 

Ablowitz,  M.J.  &:  Segur,  H.  1979  J.  Fluid  Mechanics  92,  691. 

Arkadiev,  V.A.,  Pogrebkov,  A.K.  &  Polivanov,  M.C.  1989  Physica  D36,  189. 

Boiti,  M.,  Leon,  J.J.-P.,  Martina,  L.  &  Pempinelli,  F.  1988,  Physics  Letters 
A  132,  432. 

Boiti,  M.,  Leon,  J.J.-P.  &  Pempinelli,  F.  1989,  Physics  Letters  A  141,  101. 

Boiti,  M.,  Pempinelli,  F.  &  Sabatier,  P.C.  1992  First  and  Second  Order 
Nonlinear  Evolution  Equations  from  an  Inverse  Scattering  Problem, 
(preprint). 

Borisov,  A.  et  al.  1992a  Phys.  Rev.  Lett.  65,  1753. 

Borisov,  A.  et  aJ.  1992b  Phys.  Rev.  A  45,  5830. 

Djordjevic,  V.  &  Redekopp,  L.  1977  J.  Fluid  Mechanics  79,  703. 

Fokas,  A.S.  &  Ablowitz,  M.J.  1983  Phys.  Rev.  Letters  51,  7. 

Fokas,  A.S.  &  Ablowitz,  M.J.  1984  J.  Math.  Phys.  25,  2494. 

Ichikawa,  Y.,  Mitxihashi,  T.,  Kornio,  K.  1976  J.  Phys.  Soc.  Japan  41,  1382. 

Karpman,  V.  1990  “Nonlinear  Interaction  between  Short  and  Long  Waves,” 
in  Compte  &  Boccara,  “Partially  Integrable  Evolution  Equations  in 
Physics,”  83. 

Karpman,  V.  &  Washimi,  H.  1977  J.  Plasma  Physics  18,  173. 

Kates,  R.  &  Kaup,  D.  1988  Astron.  Astrophys.  206,  9. 

Kates,  R.  &  Kaup,  D.  1989  J.  Plasma  Physics  42,  507  (Paper  I). 

Kates,  R.  &  Kaup,  D.  1989  J.  Plasma  Physics  42,  521  (Paper  II). 

Kates,  R.  &  Kaup,  D.  1991  J.  Plasma  Physics  46,  85  (Paper  III). 

Kates,  R.  &  Kaup,  D.  1992a  J.  Plasma  Physics  48,  119  (Paper  IV). 

Kates,  R.  &  Kaup,  D.  1992b  J.  Plasma  Physics  48,  397  (Paper  V). 

Kates,  R.  &  Kaup,  D.  1993,  to  be  published. 

Kaup,  D.J.  1980  Physica  ID,  45. 

Kelley,  P.  1965  Phys.  Rev.  Lett.  15,  1005. 

Shtikla,  P.  &  Stenfio,  L.  1984  Phys.  Rev.  A30,  2140. 

Spatschek,  K  1977  J.  Plasma  Physics  18,  293. 

Taniuti,  T.  1974  Suppl.  Progr.  Theor.  Phys.  55. 


20 


Zakhau'ov,  V.E.  1980,  Article  in  "Solitons”,  R.K.  BuUough  and  P.J,  Caudrey 
(eds)  [Springer- Verlag,  NY],  (see  page  250). 

Zakhau’ov,  V.E.  1972  Soviet  Physics  JETP  35,  908  (Zh.  Eksp.  Teor.  Fiz.  62 
1745  ). 

Zakharov,  V.E.  Sc  Synakh,  V.S.  1976  Soviet  Physics  JETP  41,  465  (Zh.  Eksp. 
Teor.  Fiz.  68  940  ). 

Zakharov,  V.E.  Sc  Schulman,  E.I.  1990,  Article  in  "What  is  Integrability?” , 
V.E.  Zakharov  (ed.)  [Springer- Verlag,  NY],  (see  page  230). 


21 


